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ABSTRACT 

We present an integral field spectroscopic study of the central 2x2 kpc 2 of 
the blue compact dwarf galaxy Mrk 409, observed with the Potsdam MultiAper- 
ture Spectrophotometer. This study focuses on the morphology, two-dimensional 
chemical abundance pattern, excitation properties and kinematics of the ionized 
interstellar medium in the starburst component. We also investigate the na- 
ture of the extended ring of ionized gas emission surrounding the bright nuclear 
starburst region of Mrk 409. PMAS spectra of selected regions along the ring, in- 
terpreted with evolutionary and population synthesis models, indicate that their 
ionized emission is mainly due to a young stellar population with a total mass of 
~ 1.5 x 10 6 Mq, which started forming almost coevally ~ 10 Myr ago. This stellar 
component is likely confined to the collisional interface of a spherically expand- 
ing, starburst- driven super-bubble with denser, swept-up ambient gas, ~ 600 pc 
away from the central starburst nucleus. The spectroscopic properties of the 
latter imply a large extinction (Cup > 0.9), and the presence of an additional 
non-thermal ionization source, most likely a low-luminosity Active Galactic Nu- 
cleus. Mrk 409 shows a relatively large oxygen abundance (12 + log(0/H) ~ 8.4) 
and no chemical abundance gradients out to R ~ 600 pc. The ionized gas kine- 
matics displays an overall regular rotation on a northwest-southwest axis, with 
a maximum velocity of 60 km s -1 ; the total mass inside the star-forming ring is 
about 1.4 x 1O 9 M . 

Subject headings: galaxies: dwarf galaxies - galaxies: starburst - galaxies: com- 
pact - galaxies: stellar populations - galaxies: individual (Mrk 409) 



1 Based on observations obtained at the German-Spanish Astronomical Center, Calar Alto, operated 
by the Max-Planck-Institut fur Astronomie Heidelberg jointly with the Spanish National Commission for 
Astronomy. 
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1. Introduction 

Blu e Compact Dwarf (BCD ) galaxies are systems undergoing violent bursts of star for- 



mation fjSargent fc Searldll970l ). that appear compact in the optical (starburst diameter < 1 
kpc), and have low intrinsic luminosities (Mb > — 18 mag) and low gas-phase metallicities 
(7.1 < 12 + log (O/H) < 8.3). These characteristics make them excellent nearby laboratories 
for investigating some of the most outstanding questions in the contemporary extragalactic 
astronomy: i) they are suitable sites to study the process of galaxy formation and evolution, 
as they are thought to be the loc al analogs of the buildi ng; blocks from which larger sys- 



tems were formed at high redshift (IKauffmann et al.lll993l ); ii) they allow the determination 



fPaeel et al. 


1992; 


Izotov et al. 


2007) 



galaxies, BCDs neither can sustain spiral density waves, nor suffer from disk instabilities, 
hence they permit the investigation of the star- format ion process in a relatively simple en- 
vironment. 

In order to properly characterize the BCD galaxy class, and to get insights into the above 
topics, detailed spectrophotometric analysis of individual objects, aimed to disentangle their 



ative ( 


see, e.g., 


Cairos et al. 


2002 




200' 


7). However, very few such studies have been done 


so far 


(Papaderos et al. 


2002 


2006; 


Noeske e 


t al. 


2000; 


Gil de Paz et al. 


2000; 


Fricke et al. 


2001 


: ( 


3useva et al. 


|200 


ICairos et al. 


2002|, 


2001 


). These works, all based on traditional 



observing techniques (broad- and narrow-band imaging plus long-slit spectroscopy), have 
shown that combining spectroscopy and surface photometry is indeed essential for tackling 
the problem of BCD evolution. However, they suffer from several drawbacks, the most se- 
vere being the large amount of observing time that they require: acquiring images in several 
broad- and narrow-band filters, in addition to a series of long-slit exposures sweeping the 
region of interest, translate into prohibitively long observing times of two or more nights per 
galaxy. This usually means that the data are taken under varying atmospheric conditions 
and, potentially, slight differences in the instrumental setup. Clearly, uncertainties in the 
combination of such a data set will propagate throughout in the data analysis and inter- 
pretation. Long-slit spectroscopy has also the additional problem of the uncertainty on the 
exact location of the slit. 

The relatively new observational method of Integral Field Spectroscopy (IFS) over- 
comes these probl ems. Detailed IFS studies have been carried out for SBS 0335-052 E with 
VLT /GIRAF FE Jlzotov et al.l bood) . for the low-luminosity BCD II Zw 70 wit h CAHA 



3.5m /PMAS (IKehrig et al.ll2008l ). for fiv e luminous BCDs w ith WHT/INTEGRAL (IGarcia-Lorenzo et al 
20081 ). for UM 408 with GMOS-IFU Lagos et all I2OO9I ) and for Mrk 1418 with PMAS 
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( ICairos et al.ll2009l ). Although the scope of these studies is somehow limited by a small field 
of view (FOV) or by low spectral resolution, they convincingly demonstrate that IFS is the 
method of choice for small, yet complex star-forming (SF) objects like BCDs. Each single 
exposure of an integral field spectrograph provides both spatial and spectral information, 
making the observations an order of magnitude more efficient than any traditional observ- 
ing technique. Furthermore, IFS provides simultaneously spectra for all spatial resolution 
elements, under the same instrumental and atmospheric conditions, which guarantees the 
homogeneity of the dataset. 

Systems with a high degree of asymmetry in their SF component, such as BCDs, hosting 
two or more well-separated SF regions, can be studied in a comprehensive manner only 
with 2D spectroscopy. Long standing questions in the field, as the propagation of the star- 
formation, the generation and dispersal of heavy elements or the abundance gradients, can 
be approached efficiently only with IFS instruments. 

While existing works indicate that the SF activity in BCDs proc eeds largely in bursts 



^e.g., iKriiger &: Fritze-v. Alvenslebenl Il994j ; iMas- Hesse fc Kunthlll999l ). we are still lacking 
an understanding of the two-dimensional star- formation pattern. Established is solely that, 
in the majority (90%) of BCDs, the SF regions are not randomly scattered within the low 
surface bright ness (LSB) host, but rather confined to its central part, out to ~ 2 exponential 
scale lengths (IPapaderos et al.lll996l ; ICairos et al.ll200ll ). This fact suggests a physical link 
between the mode of the SF activity in BCDs and the shape of the gravitational potential 
formed by the old, higher mass-to-light (M/L) ratio stellar LSB host. 

Some of the most metal-p oor BCDs, though t to be still in the proc ess of their forma 
tion, such as SBS 0335-052 E (F ~ 



( jThuan et al.lll999l ; iGuseva et al 



man et al.l 



2003; see 



1997 



Papaderos et al.l ll998h or SBS 1415+437 



Papaderos et al.ll2008l for a review), present clear 



signatures of uni-directional SF propagation. This process can lead to a cometary morphol- 
ogy, with the "comets head" (the dominant SF region) located at the tip of an elong ated 
LSB host delineating the trail of past propagating SF activities ( IPapaderos et al.ll2008l ). An 
irregular SF propagation pattern in form of a random-walk was reported by a detailed ana l- 
ysis of individual stellar clusters in the luminous BCD ESO 338-IG04 (jOstlin et al.ll2003l ). 
while the po t entia l role of SF propa g ation in Mrk 370 and Mrk 35 has been discussed in 
Cairos et al.l (120021 ) and lCairos et al.l (120071 ). respectively. 



In the two BCDs Mrk 86 and Mrk 409, the ring- morphology of the SF regions suggests 
that the SF activity could have been triggered by the shock wa ye produced by multiple SNs 
in a central starburst. In the case of Mrk 86, extensive studies ( )Gil de Paz et al.ll2000l . 120021 ) 
support this hypothesis. However, be cause of the smal l spatia l and spectral coverage of the 
available spectroscopic observations (jGil de Paz et al.l l2003al . hereafterGdP03), conclusive 
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results for Mrk 409 can not be drawn. 

In this paper we show the great potential of IFS when applied to BCDs by presenting 
a pilot study focused on the galaxy Mrk 409 (= NGC 3011). Mrk 409 is a lu minous (M B = 



— 17.73) BCD, belonging to the nE morphological class (lLoose &: Thuanlll986l ). Its intriguing 
gas ionized morphology — a compact nuclear starburst region surrounded by two SF rings 
at projected galactocentric radii of 5" and 18", corresponding to 0.64 and 2.3 kpc — makes it 
an excellent laboratory for investigating the star-formation pattern and abundance gradients 
in BCDs. We observed the galaxy with the PMAS spectrograph mounted on the 3.5m Calar 
Alto telescope. The PMAS field of view (2x2 kpc 2 at the Mrk 409 distance) allows us 
to map the nuclear region and the inner SF ring. The basic parameters of the galaxy are 
summarized in Table [TJ 

The paper is organized as follows: in Sect. [2] we describe the observations, the data 
analysis and the method employed to derive the two dimensional maps. In Sect. [3] we 
present continuum and emission line intensity, line ratio and velocity maps as well as results 
from the analysis of the integrated spectra of the SF regions identified in the galaxy. The 
results are discussed and summarized in Sect. [Hand Sect. 0, respectively. 



2. The Data 

2.1. Observations and Data Reduction 

Mrk 409 was observed in 2007 March, with the Potsdam MultiAperture Spectropho- 
tometer (PMAS) attached to the 3.5m telescope of the Observatorio Astronomico Hispano 
Aleman Calar Alto (CAHA). PMAS in an integral field spectrograph, with a lens array of 
16" x 16" (2x2 kpc 2 at the Mrk 409 distance) square elements. This array is connected to 
a bundle of fiber optics, whose 256 fibers are re-arranged to form a pseudo-slit in the focal 
plane of the fiber spectrograph. In the configuration used each fiber c orresponds to a spatial 



sampling on the sky of 1" x 1" (see iRoth et al.ll2005l ; iKelz et al.ll2006l for more details about 



the instrument). A grating with 300 grooves per mm was used during the observations, in 
combination with a 2048 x 4096 pixel SITe ST002A CCD detector; this provides a spectral 
range of 3590-6996 A with a linear dispersion of 3.2 A per pix and a resolution (FWHM) 
of about 7 A. 

We observed a total of 6000 sec on the galaxy, with the integration time split into five 
exposures of 1200 sec each; additional sky frames were taken moving the telescope several 
arcmin away from the target position. Calibration frames were observed before and after the 
galaxy exposures. These calibrations consist of spectra of emission line lamps (HgNe lamp), 
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Table 1. Basic data for Mrk 409 



Parameter 


Data 


Note 


Other names 


NGC 3011, UGC 5259 




RA (2000) 


09 49 41 




DEC (2000) 


+32 13 16 




Wheiio (km s" 1 ) 


1527 




U (MpcJ 


ZO.O 




A B (mag) 


0.071 


(1) 


M B (mag) 


-17.73 




rriB (mag) 


14.37 ±0.03 


(2) 


m R (mag) 


13.33 ±0.05 


(2) 


F(B.a) (ergs s -1 cm" 2 ) 


2.93 ±0.17 x 10~ 13 


(2) 


mj (mag) 


11.88 ±0.02 


(3) 


van (mag) 


11.34 ±0.03 


(3) 


m Ks (mag) 


11.17 ±0.04 


(3) 


M HI (M ) 


0.31 x 10 9 


(4) 


M T (M ) 


0.33 x 10 10 


(4) 



Note. — Coordinates, heliocentric velocity and distance, 



all taken from NED ( |http:/ /nedwww.ipac.caltech.edu/ ). 
Distance calculated using a Hubble constant of 73 km s _1 
Mpc -1 , and taking into account the influence of the Virgo 
Cluster, the Great Attractor and the Shapley superclus- 
ter. Absolute magnitude in the B band, computed from 
the tabulated integrated B magnitude and distance. (1) 
Absor ption coefficient in the B band , from ISchlegel et al. 
Jl998h . (2) iGil de Paz etall J2003bh . (3) From 2 MASS. 
(4) Neutral hydrogen m ass Mhi and total mass Mr from 
Thuan fc Martini (jl98lh . 
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which are required to perform the wavelength calibration, and spectra of a continuum lamp, 
necessary to locate the fiber spectra on the CCD and to perform the flat-fielding correction. 
The spectrophotometric standard stars BD+75325 and BD+332642 were observed for flux 
calibration. As usual, bias and sky-flat exposures were taken at the beginning and at the 
end of the night. The seeing was about 1.8 arcsec. 

The data were processed using standard IRAF@ tasks. The reduction procedure included 
bias subtraction and image trimming, tracing and extraction, wavelength and distortion cal- 
ibration, flat-fielding, combination of the individual frames, sky subtraction and flux calibra- 
tion. For a d e tailed description of the individual steps performed in fiber data processing see 
Cairos et al.l ( 20091 ). The typical errors on the wavelength calibration are 0.02 A. The flux 
calibration was carried out using our observations of spectrophotometric standard stars. By 
comparing the sensitivity curves obtained for all the spectrophotometric standards observed 
throughout the three nights of the observing run, we estimate that the relative uncertainty 
on the calibration factor is generally equal or less than 2%, except blueward of 4000 A, where 
the curves show a marked change of slope and the uncertainty increases up to about 8%. 



2.2. Line Fitting and Map Generation 

In this study we generate 2D intensity maps for the most important emission lines as well 
as for different continuum regions. The relevant parameters of the emission lines (center, flux 
and equivalent width) were determined by fitting single Gaussians; the fit was carried out 
by using the y 2 minimization algorithm implemented by C. B. Markwardt in the mpfitexpr 



IDL library (IMarkwardtl 120091 ) . The H/3 line, in which the absorption component was also 
visible, was modeled as the sum of two Gaussians. The continuum (typically 30-50 A on 
both sides) was fitted by a straight line. Lines in a doublet were fitted imposing that they 
have the same redshift and width. 

Criteria such as flux, error on flux, velocity and width were used to do a first, automatic 
assessment of whether to accept or reject a fit. For instance, lines with too small (less than 
the instrumental width) or too large (> 5 A) widths were flagged out, similar to lines with 
an error on the flux of more than about 10% (the exact limits depending on the specific 
line). In order to improve the reliability and stability of fits invoking both an absorption 
and emission component, an additional criterion was applied to the equivalent width of the 



2 IRAF is distributed by the National Optical Astronomy Observatories, which are operated by the As- 
sociation of Universities for Research in Astronomy, Inc., under cooperative agreement with the National 
Science Foundation. 
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absorption component (allowed to vary between and 5 A). 

Each individual fit was visually checke d, and when necess ary, the aforementioned auto- 
mated quality criteria were overridden (see lCairos et al.ll2009l ). 

The fit procedure gives, for each spectral line, a table with the spaxel number and 
the values of the measured parameter; this table is then used to build a 2D spatial map. 
Continuum maps were built in the same way by summing the signal within wavelength 
intervals free from emission lines or strong sky-lines residuals. 



3. Results 

3.1. Continuum and Emission Line Intensity Maps 

The main advantage of IFS is that it maps an object simultaneously in different contin- 
uum bands and emission lines, so that we can obtain at the same time the spatial distribution 
of the stellar component and of the ionized gas emission. We constructed continuum and 
flux maps of the brightest emission lines in Mrk 4 09. To allow a quick com parison with our 



results, the 5-band and Ha images of the galaxy fjGil de Paz et al.ll2003bl ) are displayed in 
Figure HJ 

Figure |2j displays the [O II] A3727, H/5, [O III] A5007, [O I] A6300, Ha, [N II] A6584 
and [S II] AA6717, 6731 emission-line maps. These maps trace the regions of recent star 
formation. In all maps except one Mrk 409 shows a similar, irregular pattern, with a number 
of SF knots distributed in a ring of diameter ~ 5 arcsec (600 pc) around the nuclear starburst 
(region #3); the ring is resolved into at least five smaller SF knots. Only in [O I] A6300 
the galaxy has a strikingly different morphology, revealing solely a single prominent central 
peak. The presence of strong [O I] A6300 is an indicator that power-law photoionization or 
shock-heating are important in the nuclear region of Mrk 409. 

The stellar continuum emission in the spectral range 4500-4700 A is also shown in 
Figure [2J This map displays a regular morphology with elliptical isophotes with nearly 
constant position angle, and a well defined central peak, whose location is the same in all 
continuum bands and also coincides with the nuclear SF region #3. 



Fig. 1. — left panel: B-band image of Mrk 409; the central box indicates the field of view 
covered by our PMAS observations; right panel: Continuum-subtracted Ha image with the 
individual SF knots labeled. The field of view is one arcmin in both fra mes. North is up, 



east to the left. Both images, retrieved from the NED, are published in iGil de Paz et al. 
f |2003bh 
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Fig. 2. — Emission line intensity maps derived from Gaussian fits to each line and spaxel 
(see text for details): [O II] A3727; H/J; [O III] A5007; [O I] A6300; Ha; [O III] A5007; 
[S II] A6717, 6731. The last panel shows the continuum emission between 4500 and 4700 A, 
a spectral region free from line emission ("pure continuum"). North is up, east to the left. 
Axis units are arcseconds. All images are shown in logarithmic scale to bring up the fainter 
regions. Flux units are 10~ 18 ergs cm~ 2 s _1 A" 1 . 
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3.2. Line Ratio Maps 

In order to investigate the main ionization mechanisms operating in Mrk 409, we built 
the [O III] A5007/H/3, [O I] A6300/Ha, [N II] A6584/Ha and [S II] AA6717, 6731/Ha line 
ratio maps, shown in Figure [3j In these maps high excitation corresponds to high values 
of [O III] A5007/H/3, and to low values of [N II] A6584/Ha and [S II] AA6717, 6731/Ha. 
All the excitation maps display a similar pattern. In the SF knots in the ring, the values 
of the different ratios are those characteristic of H II regions. In the central starburst, 
however, opposite to what is expected in regions photoionized by stars, high values for 
all ratios are found: [O III] A5007/H/3> 6; [O I] A6300/Ha > 0.2; [N II] A6584/Ha > 
0.5; [S II] AA6717, 6731/Ha > 1), most probably indicating the presence of an additional 
ionization mechanism — an Active Galactic Nucleus (AGN) or shocks. 

The extinction in the galaxy (Figure[3]) has been derived from the Ha/H/3 ratio. Mrk 409 
displays an inhomogeneous extinction pattern, which peaks in the nuclear starburst. In this 
region the Ha/H/3 line ratio rises to values as high as 6, suggesting strong extinction by dust 
with a reddening E(B — V) of up to ~ 0.70 mag. By contrast, in the c ircumnuclear SF knots 



the H a/H/3 ratio is only slightly larger than the nominal value of 2.86 (jOsterbrock fc Ferland 



20061 ). and implies E(B — V) < 0.2. Note that, applying the extinction derived in the 
central region to the whole SF component, as is usually done in long-slit or single-aperture 
spectroscopic studies (with, e.g., the SDSS), would lead to a large overestimate of the Ha 
luminosity and consequently of the Star Formation Rate (SFR). In this respect, a spatially 
resolved 2D determination of the extinction is essential for a proper characterization of the 
stellar populations and the star formation history of Mrk 409, and more in general of BCD 
galaxies. 

In computing the extinction map, we corrected the H/3 line for underlying stellar ab- 
sorption by fitting a Gaussian profile to its absorption wings. For the Ha line, the absence 
of absorption wings makes a reliable decomposition impossible. In this case, in principle 
several approaches could be adopted : a) to take for Ha the same value in absorption that 
for H/9, that is, to set EW(Ha a b s ) = EW(H/3 abs ), traditionally the most common strategy; 
b) to apply spectral synthesis models to disentangle the young and old population spectra 
and, from the model output, estimate the EW(Ha) in absorption; c) to assume that the 
absorption in Ha is negligible. 

Taking into account that both a) and b) rest on "model-dependent" values, probably 
overestimations of the underlying absorption in Ha, which in turn overestimates Ch/3, we pre- 
fer to adopt the most conservative c) approach, keeping in mind that under this assumption 
the computed extinction value is a lower limit to the actual value. 
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Fig. 3. — Ionization ratios (first four panels), interstellar extinction and electron density 
maps. Axis units are arcseconds; north is up, east to the left. Ionization ratio maps are in 
logarithmic scale. 
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3.3. Integrated spectroscopy of selected galaxy regions 

As the next step in our analysis, we used the reconstructed flux and continuum maps 
to delimit the regions of interest in the galaxy. A total of six regions were identified, the five 
SF located in the ring and the nuclear SF region #3 (see Fig. 

We identified the spaxels corresponding to each of these regions, and summed their 
spectra together to create the integrated spectrum. In this way we produce higher signal-to- 
noise spectra (in comparison to the spectra of the individual spaxels), better suited for the 
derivation of physical parameters and chemical abundances. 

In Figure H] we compare the flux-calibrated spectra of these regions, prior to correction 
for interstellar extinction. The spectra are arbitrarily shifted on the vertical axis, for the 
sake of clarity. It can be seen that the spectral energy distribution (SED) of the nucleus 
#3 differs strikingly from that in regions along the circumnuclear SF ring (#1, #2, #4-6). 
In the nuclear starburst the stellar continuum is an order of magnitude higher than in the 
ring; it displays prominent absorption features, including the higher-order Balmer lines, H<5 
and H7, Call H and K lines, Nal A5922 or the G-band at A4304. Its strong [O I] A6300 
emission is also remarkable, suggesting an additional contribution to the ionization, other 
than photoionization from OB stars. This comparatively red stellar continuum and the strong 
absorption features both indicate that the emission is dominated by older stars. The SF 
regions in the ring display, on contrast, an almost flat continuum with stronger emission lines, 
such as [O II] A3727, H/3, [O III] A5007, Ha, [N II] AA6548, 6584 and [S II] AA6717, 6731, 
in addition to weak [O I] A6300 emission. These spectra also reveal absorption wings in H<5 
and H7 and such absorption features as Call H K AA3933, 3969. 

We measured the fluxes and equivalent widths of the most prominent emission lines in 
each spectrum. The extinction in each region was determined from the Ha/H/3 ratio, as 
the measured H<5 and H7 fluxes are doubtful due to their intrinsic weakness and the large 
uncertainties in the correction for the underlying stellar absorption. For the theoretical 
Ha/H/3 value, we adopted the ca se B Balmer recom binat ion decrement H a/H/3 = 2.86 for 



T = 10000 K and N P = 10 4 cm" 3 (iBrocklehurstl 119711) : the IWhitfordl (119581 ) reddening curve 



as parametrized for iMiller fc Mathewsl (119721 ) was used 



For H/5, where absorption wings are clearly resolved, we fitted simultaneously an ab- 
sorption and an emission component, in order to obtain precise measurements of the H/3 flux 
in emission; however, in Ha, the absence of absorption wings makes a reliable decomposition 
impossible. Again, we assume that the absorption in Ha is negligible, so that the derived 
Ch/? are actually lower limits. 

For comparison, we also tabulate the values obtained by setting EW(Ha a b s ) = 2 A, as 
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Fig. 4. — Flux-calibrated spectra of the six resolved star-forming regions in Mrk 409. The 
spectra are shown in logarithmic scale and are offset for the sake of clarity. 
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in the output spectra from the synthesis models run in Sect. [U we measured an EW(Ha a b s ) 
of about 1.5-2 A. For this latter case we see a modest increase in the extinction coefficient. 

The reddening corrected line intensities in each region are listed in Table [2J The most 
relevant line ratios are shown in Table [31 

Th Ha flux emitted by the SF regions is about 60% of the total Ha flux within the 
PMAS field of view. 



3.3.1. Constraining the Ionization Mechanism 



The ionization mechanisms acting in a galaxy and their spatial distri bution can be 



studied by means of the line rati o diagrams for the usual diagnostic lines ( iBaldwin et al. 



198ll ; IVeilleux fc Osterbrockl 119871 ) . These diagnostic diagrams help distinguish among the 
different ionization mechanisms (photoionization from young stars, photoionization by a 
power law continuum source and shock-wave heating) operating in a galaxy. 

The position of the six regions under study on the diagnostic diagram [O III] A5007/H/3 
vs [O I] A6300/Ha, [N II] A6584/Ha and [S II] AA6717 6731/Ha is shown in Figur e El The 
empirical boundaries between the different zones (fromlVeilleux fc Osterbrock!ll987l ). as well 
as the theoretical boundaries proposed by iKewley et al.l (120011 ). are also plotted. 



H II-like ionization (i.e. star formation) seems to be the dominant mechanism in the 
five SF knots conforming the ring (see also Table [3] for the specific values of the emission 
line ratios). The line ratios of the nuclear region #3, however, are consistent with those of 
Seyfert galaxies. We came to the same conclusion assuming an extreme correction for the 
underlying Ha absorption of EW(Ha a b s ) = 4 A for the nuclear region. The possibility of 
the nuclear region of Mrk 409 containing a non-thermal energy source cannot be discarded. 
Indirect evidence for this interpretation comes from the failure of purely stellar synthetic 
spectra to reproduce the observed spectrum of that region (see Sect. H]). 



3.3.2. Physical parameters and chemical abundances 

The physical properties and oxygen abundances for Mrk 409 are listed in Table [31 

The electron density (N e ) for the individual SF knots was derive d from the [S II] AA6717, 6731 
A line ratio using the task temden in the IRAF nebular package (IDe Robertis et al.lll987l ; 



Shaw &: Dufourlll995[ ). 



Table 2. Reddening corrected line intensity ratios 



Line 


Ion 




Knot #1 


Knot #2 


Knot #3 


Knot #4 




Knot #5 






Knot #6 


( A ) 






Fx 


-Wx 


F X 


-Wx 


Fx 


-Wx 


F X 




-Wx 


Fx 




-Wx 




Fx 


-Wx 


3727 


[Oil] 


5 


06 ± 1.18 


40.3 ±2.9 


4.29 ±0.80 


46.8 ± 2.6 


13.13 ±2.61 


25.5 ±0.7 


4.84 ± 0.62 


52. 


3 ±1.8 


5.03 ± 0.68 


47 


,0± 1.9 


3. 


50 ±0.34 


86.1 ± 1.2 


4861 


H/3 


1. 


00 ±0.00 


5.4 ±0.7 


1.00 ±0.00 


7.9 ±0.8 


1.00 ±0.00 


1.4 ±0.2 


1.00 ±0.00 


9 


1±0.5 


1.00 ±0.00 


7 


,4 ±0.5 


1. 


00 ± 0.00 


18.6 ±0.3 


4959 


[o HI] 





89 ±0.15 


5.2 ±0.4 


0.48 ± 0.07 


3.9 ± 0.4 


2.51 ±0.34 


3.8 ±0.1 


0.61 ±0.05 


5 


9 ±0.3 


0.74 ±0.07 


6 


,0±0.4 


0. 


79 ± 0.03 


14.3 ±0.2 


5007 


[o III] 


2. 


13 ±0.30 


12.6 ±0.4 


1.64 ±0.18 


13.7 ±0.4 


5.68 ±0.76 


8.9 ±0.1 


1.89 ±0.12 


18. 


6 ±0.4 


1.98 ± 0.14 


16 


.3 ±0.4 


2. 


36 ±0.07 


43.4 ± 0.4 


6300 


[O I] 





20 ± 0.06 


1.3 ±0.3 


0.16 ±0.03 


1.5 ± 0.2 


0.44 ± 0.09 


1.0 ±0.1 


0.18 ±0.03 


2 


1±0.3 


0.17 ±0.03 


1 


,7 ±0.3 


0. 


10 ±0.01 


2.0 ±0.3 


6548 


[N II] 





25 ±0.07 


1.7±0.3 


0.25 ±0.05 


2.5 ± 0.3 


0.40 ± 0.09 


1.0 ±0.1 


0.27 ±0.03 


3 


.3 ±0.3 


0.25 ± 0.03 


2 


.6 ±0.2 





21 ±0.01 


4.5 ±0.2 


6563 


Ha 


2 


86 ±0.55 


19.7 ±0.4 


2.86 ± 0.43 


28.9 ± 0.6 


2.86 ±0.52 


7.5 ±0.1 


2.86 ±0.24 


35 


6 ±0.6 


2.86 ± 0.27 


29 


,9 ±0.4 


2 


86 ±0.12 


60.1 ±0.3 


6584 


[N II] 





58 ±0.13 


4.0 ±0.3 


0.62 ±0.10 


6.2 ± 0.4 


1.19 ±0.22 


3.1 ±0.1 


0.67 ±0.06 


8, 


3 ±0.3 


0.62 ± 0.06 


6 


,4 ±0.3 


0. 


53 ±0.02 


6.4 ±0.3 


6717 


[S II] 





58 ±0.13 


4.2 ±0.4 


0.49 ± 0.08 


5.1 ± 0.3 


0.89 ±0.17 


2.5 ±0.1 


0.67 ±0.06 


8, 


8 ±0.3 


0.57 ±0.06 


6 


,2 ±0.4 


0. 


40 ±0.02 


0.0 ±0.4 


6731 


[S II] 





44 ±0.10 


3.2 ±0.3 


0.35 ±0.06 


3.7 ±0.3 


0.86 ±0.16 


2.4 ±0.1 


0.51 ±0.05 


6 


7 ±0.4 


0.39 ± 0.05 


1 


,3 ±0.3 


0. 


30 ±0.02 


0.0 ±0.3 



Cup 0.23 ±0.19 0.22 ±0.14 0.91 ±0.18 0.21 ±0.08 0.27 ±0.09 0.21 ±0.04 

Cup* 0.35 ±0.19 0.30 ±0.14 1.23 ±0.18 0.27 ±0.08 0.35 ± 0.09 0.22 ±0.04 

F(K(3) 1.74 ±0.84 2.51 ± 0.94 18.27 ± 8.91 4.09 ± 0.86 3.63 ± 0.87 117.27 ±1.74 



as 

Note. — Reddening-corrected line intensities, normalized to F(H/3) = 1, for the individual regions selected in Mrk 409. The reddening coefficient, Ch/3, and the corrected H/3 flux, F(H/3) 
(xl0 — 15 ergs cm -2 s _1 ), were computed assuming no underlying Ho absorption. For comparison, we also list Ch/3*, the reddening coefficient derived by assuming EW(Ha a j, s ) =2 A. 
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Fig. 5.- 



LINERs from 


H 


II 


Veilleux & Osterbrock 


Kewlev et al. 


( 


2001 


)• 



Optical emission-line diagnostic diagrams, separating Seyfert galaxies and 
igion-like objects. Solid lines are the empirical borders from 
(119871 ). while dotted lines are the theoretical borders from 



Kewlev et all (120011 ). The line ratios for the nuclear SF region of Mrk 409 (#3) and the 



five circumnuclear SF regions (#1-2 and #4-6) are plotted with open and filled symbols, 
respectively. For the nuclear region, we also mark (open star without errorbar) the position 
corresponding to assuming a correction for underlying Ha absorption of EW(Ha a b s ) = 4 A 
(for all other knots, this shift would be unnoticeable) . Note that, in all three classification 
diagrams, the nuclear SF region falls in the Seyfert regime. 
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The " direct" oxygen abundance method relies on measurements of temperature-sensitive 
auroral lines (usually [O III] A4363), which are weak and often not detected in oxygen-rich 
systems (above ~ 12 + log(0/H) = 8.3 — 8.4). This is the case of the spectra of Mrk 409, in 
which either we fail to detect the [O III] A4363 emission line, or, in the cases where we do 
detect the [O III] A4363 line, the S/N is too low for a reliable flux determination. 

Different empirical and theoretical metallicity calibrations have been devised over the 



past two decades as an alternative 



or estimating oxygen abundances without any direct 



electron temperature measurements ( 


Paeel et al. 


19791: 


Alloin et al. 


19 


79|; 


Edmunds &; Paeel 


1984; 


McGaueh 


1991; 


Vilchez & Esteban 


1996: 


3 ettini & Paeel 


2004; 


Pilvuein et al. 


2006). 



These calibrations are referred to as strong-line methods since they are based on ratios built 
on strong emission lines ([O II] A3727, [O III] AA4959, 5007, [N II] A6584JS I I I AA6 717, 6731, 
[S III] AA9 069, 9532). We have us ed the strong-line methods from lMcGauehl (1 1 99 ll ) [hereafter 
M91] and iPettini fc Paeell (120041 ) [hereafter PP04] to derive the oxygen abundances of all 
selected SF knots, except the nuclear region #3 (see below). 

Efforts have been made to refine the calibration of the R23 parameter. One of the most 
successful is the calibration by M91, which is based on pho toionization models — the analytic 
expres sions for the M91 theoretical models can be found in Kobulnickv. Kennicutt fc Pizaeno 
(119991) . This calibration can be as precise as 0.2 dex (IKobulnicky. Kennicutt fc Pizaeno 
1999h . However, due to the well-known double- valued behavior between R23 and O/H, this 
method requires some prior knowledge of the galaxy's metallicity, in order to place it on the 
correct branch of the curve. We used the ratio [N II] A6584/Ha to break the R23 degeneracy; 
in all the SF regions, log([N II] A6584/Ha) is > —1.1, indicating that the knots lie on the 
upper branch — the division between the R23 upper and lower branches takes place between 
-1.3 < log([N II] A6584/Ha) < -1.1 fcewlev fe Ellison! boosh . 



On the other hand, PP04 revised the N2 (N2 = log[N II] A6584/Ha) and 03N2 (03N2 
= log([0 III] A5007/H/5/[N II] A6584/Ha)) indices, using 137 extragalactic HII regions. The 
uncertainties in metallicity determinations based on these calibrations amount to ~ 0.38 and 
0.25 dex, respectively. 

Inspection of Fig. [5] shows knot #3 to be located in the Seyfert zone of the diagnostic 
diagram. For this region, we estimate therefore the metal licity using the chemical abundance 
calibrations proposed by IStorchi-Bergmann et al.l (119981 ). These calibrations were derived 
for the range 8.4 < 12 + log(0/H) < 9.4 and tested using a sample of Seyfert galaxies and 
LINERs that have H II regions in the vicinity of their nucleus. 



Table 3. Line ratios, physical parameters and chemical abundances 



Parameter 


Knot #1 


Knot #2 


Knot #3 


Knot #4 


Knot #5 


Knot #6 


log( [ ° n ^ A5007 ) 

, / [N II] A6584 \ 

, / [S II] AA6717, 6731 n 
log( l Hq ) 

log(^^oo) 


0.33 ± 0.06 
-0.70 ±0.09 
-0.45 ± 0.08 
-1.16 ± 0.13 


0.22 ± 0.05 
-0.66 ± 0.07 
-0.53 ± 0.06 
-1.26 ± 0.10 


0.75 ± 0.06 
-0.38 ±0.08 
-0.21 ±0.07 
-0.81 ±0.09 


0.28 ± 0.03 
-0.63 ±0.04 
-0.38 ±0.04 
-1.21 ±0.07 


0.30 ± 0.03 
-0.67 ±0.05 
-0.48 ± 0.04 
-1.23 ±0.09 


0.37 ± 0.02 
-0.73 ±0.03 
-0.61 ±0.02 
-1.46 ±0.06 


N c 

12 + log(0/H) a 
12 + log(0/H) 6 
12 + log(0/H) c 
12 + log(0/H) d 
log(N/0) e 


117 
8.40 
8.46 
8.40 

-1.13 


< 100 
8.45 
8.49 
8.53 

-1.02 


516 

8.44 
-1.28 


102 
8.44 
8.51 
8.45 

-1.05 


< 100 

8.48 
8.42 

-1.11 


< 100 
8.38 
8.43 
8.54 

-1.01 



Note. — The line ratios have been corrected for galactic extinction. The quoted uncertainties include both measurement 
errors and the uncertainty in the calibration factor, (a) O xygen abundances com puted using the fit to the relationship 
between T e metallicities and the ([O III]/H/3)/([N II]/Ho) in lPettini fc Pa gcl (2001); (b) oxygen a bundances obtained from 
the relationship b etween T e metallic ities and the [N II] /Ho ratio given bv fpettini h Pagell l|2004h : (c) oxygen abundances 
derived using the iMcGaughl <|l99lj) calibra tion for the upper branch; (d ) oxygen abundance computed from chemical 
abundance calibr ations for Seyfert galaxies l|Storchi-Bergmann et al.lll998h ; (e) [N II] / [O II] was used to obtain the N/O 
abundance ratio l|Perez-Montero fc Di'a"3l2005h . 
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3.4. Kinematics of the ionized gas 

Figure [6] shows the velocity field of the ionized gas, derived from the [O III] A5007 
and the Ha emission lines. Both velocity maps show essentially the same regular pattern, 
although in the [O III] A5007 map it appears slightly more distorted; the galaxy displays an 
overall smooth rotation pattern along a northwest-southeast axis. 

Figure [7] shows the velocity profile and the position angle of the kinematical axis, as 
determined b y a tilted-ring analysis of the Ha velocity field based on the kinemetry method 
developed by lKrajnovic et al.l (120061 ). 



The velocity rises rapidly in the inner 4 arcsecs, then flattens out at about 60 km s -1 . 
Adopting for simplicity the virial formula valid for sphe rical symmetry , we can estimate the 
mass from the relation: M(R) = 2.32 x 10 5 Rv^ IC (R) (jLequeuxill983l ). where R is in kpc, 
^circ = v obs/ sin(i) is the circular velocity, % is the inclination angle, and M is expressed in 
solar units. 

The apparent flattening of Mrk 409 is q = b/a ~ 0.80 (SDSS, Data Release 6), which, 
for an intrinsic flattening gi ntr < 0.3 gives an inclination angle i ~ 38. Taking an observed 
rotational velocity v ~ 60 km s" 1 at the position of the SF ring, we obtain an indicative 
virial mass of ~ 1.4 x 10 9 M Q interior to the SF ring (i? r i ng ~ 0.6 kpc). 
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Fig. 6. — Velocity field of the ionized gas in the central region of Mrk 409 in the [O III] and 
Ha lines. Axis units are arcseconds; north is up, east to the left. 
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Fig. 7. — Velocity curve and position angle profile of the kinematical major axis (measured 
from north through east) as determined by a tilted-ring model analysis. 

4. Discussion 



As it was already stated in the introduction, one of the main open questions in BCD 
research is how the star formation originates and propagates across these galaxies. The 
spatial distribution of the SF knots in the central part of Mrk 409 — a nuclear starburst 
surrounded by a SF ring — makes this galaxy a paradigm case of study. 

An extended circumnuclear SF zone may play an important role on BCD evolution, 
contributing to the build-up of the stellar component and to the chemical enrichment of 
the Interstellar Medium (ISM). In addition, the possibility of Mrk 409 hosting a non- 
thermal nucl ear energy sourc e (see Sect. 13.3.11) . makes this object even more intriguing. 
In this line, llzotov fc Thuanl (120081 ) reported recently on the detection of AGN candi- 
dates in four metal-poor dwarf galaxies, and c ircumnuclear SF a c tivity is known to ex- 



ist in several AGNs hosted by normal galaxies (I Wilson et al.l Il986l ; iKotilainen et al.l 12000 
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Gonzalez Delgado et al. 1 1200 lh 



GdP03 suggested that the SF ring in Mrk 409 could have been formed as the result 
of the interaction of a starburst-driven shock with the sourrounding ISM. These authors 
performed echelle spectroscopy of the central part of Mrk 409 and found a doubly-peaked 
Ha emission line profile, with its two peaks shifted by —54 ± 3 and 37 ± 3 km s _1 relative 
to the systemic velocity of the stars (presumably because of the lower spectral resolution, 
we do not see this feature in our spectra). But their data, covering only a thin slice of the 
BCD's central part, did not allow to firmly validate this scenario. 

This idea appears also viable from the viewpoint of observational and theoretical work 
that indicates that BCDs expe rience multiple bursts of star formation in their central 



part through their lifetimes (e.g. iKriiger fc Fritze-v. Alvenslebenlll994l ; iMas- Hesse fc Kunth 



19991 ). The mechanical energy released by the collective action of stellar winds from mas- 
sive stars and SNe can readily result in large-scale expansive motions i n the ISM reflected 
on e.g. Ha super-s h ells out to kpc sc ales from the starburst region (IMeurer et al.l Il992 



Marlowe et al. 1995; Cairos et al. 



20011 ) . Observations of BCDs indicate expansive veloc 



1997 




van Evmeren et al. 


2007 


) and > 200 km 


Johnson et al. 


2000; Haro 2, 


Leerand et al. 


1997) 



Meurer et al. 1992; Marlowe et al. 



199 



5; 



Martin 



19991 : 



An alternative hypothesis, which must also be taken into account, is that the ring found 
in Mrk 409 is a resonance ring (RR), similar to those observed in barred galaxies fjButa 
19881 ; iButa fc Purcelllll998l ; iKotilainen et al.ll2000l ). RRs are thoug ht to result from efficien t 
gas inflow onto the central region of a galaxy, driven by the bar (ICombes &; Gerinl Il985l ). 
and its collapse at the radius of the inner Lindblad resonance. While RRs are typically 
observed in galaxies of normal mass, they may as well exist in lower-mass systems, such as 
Mrk 409. Note that the gas-phase metallicity in the SF ring in Mrk 409 is in the ra nge of 
values determined in circumnuclear SF regions in normal galaxies (IDiaz et al.l 120071 ) . and 
also that the SF ring compares well in size (R~600 pc) with RRs in normal galaxies - the 
ave rage radius of galactic RRs is 750 pc (IButa fc Crockerlll993l ) with values as small as ~200 
pc (jComeron et al.ll2009i ). 



In order to distinguish between the two aforementioned scenarios it is important to study 
the SFH of the individual knots along the ring. In the framework of the RR hypothesis, it 
is reasonable to assume that this feature has been stable over several orbital periods (a few 
10 8 yr), forming stars at a nearly constant rate, whereas the blast wave hypothesis, on the 
other hand, requires that the SF ring as a whole was created nearly instantaneously in the 
recent past. 
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Our PMAS data, mapping the nucleus and the SF ring, permit, despite their moderate 
spectral resolution, the investigation of this issue. In order to constrain the SFH of each 
SF knot, we interpret their integrated spectrum (Section 13.31) by means of spectral synthesis 
models. To chec k the consistency of our results, we use both the population synthesis 
code starlight (ICid Fernandes et al.ll2004 l2005allbl; iGarcia-Rissmann et al.ll2005l) and the 
evolutionary spectral synthesis code pegase ( Fioc fc Rocca-Volmerange Il997 ). 



We use starlight to synthesize the observed stellar continuum in each SF region 
as due to the superposition of single-age stellar populations (S SPs) of different ages and 
metallicities. The SSP library used is based on stellar models by lBruzual fc Charlotl (120031 ) 
for a Salpeter initial mass function (IMF), five metallicities (between Z /5O and 2.5 Z Q ) and 
59 ages (between 1 Myr and 13 Gyr). Prior to modelling, the observed spectra were corrected 
for interstellar extinction (using the Ch/3 derived without stellar absorption correction in Ha), 
and spectral regions containing strong emission lines or skyline subtraction residuals were 
flagged out. 

Fig. El (top panel) displays the extinction-corrected spectra of the circumnuclear SF 
regions #1, #2 and #4-6, normalized at 4200 A (orange color) together with the best- 
fitting synthetic spectra overlaid in black. Flagged intervals in each spectrum are marked 
with shaded vertical bars. The smaller panels to the upper-right and lower-right show the 
age distribution of the SSPs selected by starlight as a function of, respectively, their 
luminosity and mass fraction in %. We find that the spectra of all circumnuclear regions 
are well fit by the superposition of a several Gyr old stellar component with a young (<50 
Myr) population which dominates the optical light. In none of the regions do starlight 
models indicate an intermediate-age stellar population bridging the gap between 0.1 and ~1 
Gyr. All solutions imply no or little {Ay < 0.14 mag) intrinsic extinction in excess to that 
inferred from the measured Ch/3- According to the fits, the formation of the underlying host 
galaxy of Mrk 409 was accomplished between ~ 1 and ~ 10 Gyr ago, and at least 3/4 of its 
stellar mass was at place ~ 6 Gyr ago. This is in good agreement with the red B — R color 
of - 1.1 mag of the host of Mrk 409 (GdP03). 

Secondly, we apply a modified version of pegase 2.0 to study whether the SFH along 
the ring is better approximated by an instantaneous or a continuous model (models inst and 
c0-c2, respectively, in Table H]). In this approach, each of the observed spectra is modeled 
as due to the superposition of the SED of an old and a young stellar population describ- 
ing, respectively, the underlying old host galaxy and the SF ring. We use the same input 
spectra as for starlight after resampling to the lower spectral resolution of pegase. We 
furthermore assume a Salpeter IMF between 0.1 and 100 M and a fixed stellar metallicity 
of Zq/5. For the old underlying population we adopt an exponential SFR with an e-folding 
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Fig. 8. — Spectra of the five selected regions in the circumnuclear star-forming ring in 
Mrk 409, corrected for extinction and normalized by the observed flux at 4200 A (orange 
color). The spectra are arranged from region #1 (top) to region ^6 (bottom). The best- 
fitting synthetic spectral energy distribution ( SEP) for the stellar continuu m, calculated with 
the population synthesis code STARLIGHT (ICid Fernandes et al.ll2005al ). are shown in the 
upper panel in black. Flagged intervals are marked with shaded vertical bars. The smaller 
upper-right and lower-right panels show the age distribution of the stellar populations in 
the synthetic SEDs as a function of, respectively, their luminosity (at 4200 A) and mass 
contribution in %. The SEDs corresponding to the best-fitting Pegase 2.0 models inst, 
cO and cl are shown in the lower panel with the black, solid-gray and dashed-gray curve, 
respectively. The blue and red curves illustrate the young and old stellar component, as 
derived from the inst model. 



-25 - 




4000 4500 5000 5500 6000 
wavelength (A) 



6500 



— 80 
p 

^ 60 

« 40 
u 

s 20 

- o 

100 

c 10 



0.1 



'0.01 



7 8 9 
log (t. [yr]) 



10 




4000 



4500 5000 5500 6000 6500 
wavelength (A) 



log (t [yr]) 




Fig. 8 - 



Continued. 
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time r = 3 Gyr, except for model c2 for which r was set to 1 Gyr. For these two e-folding 
times the observed B-R color of the host galaxy is reproduced by a present age of 10.3 Gyr 
and 5 Gyr, respectively. Models inst and cO fit both the stellar SED and the observed EWs 
of the Ha and H/3 emission lines whereas in models cl and c2 only the stellar SED was 
fitted. 

The results for each SF region are summarized in Table HI The age and mass of the 
young stellar component tsF and M* sf are given in cols. 3 and 4, respectively. Columns 5 
and 6 list the predicted EWs of the Ha and H/3 emission line and B — R color. Column 7 
gives the difference A/i# between the 5-band surface brightness of the old stellar population 
and that of the total stellar population predicted by the fit (/xb.oH — A*s,totai), and col. 8 the 
reduced x 2 - The observed EWs of the Ha and H/3 lines are also indicated for each SF region. 

The best-fitting stellar SEDs for each circumnuclear region, overlaid with the observed 
spectrum, are shown in Fig. [8] (bottom panel). The SEDs corresponding to models inst, 
cO and cl are shown with the black, solid-gray and dashed-gray curve, respectively. The 
blue and red curves illustrate the young and old stellar component, as derived from the inst 
model. 

Inspection of Table H] and Fig. [8] shows that the instantaneous model reproduces best 
both the stellar SED and the EW of Ha and H/3, implying in all cases a burst age tsF of 
6 to 10 Myr. The predicted stellar mass from this model adds up to ~ 1.5 x 10 6 . Models 
cO, assuming continuous star formation since tgp (in most cases > 2 Gyr), generally fail 
to simultaneously reproduce both the EWs and the stellar SED, with a clear tendency for 
underestimating the latter for A < 4500 A. Continuous SF models fitting the stellar SED 
only (cl and c2) predict, similar to the instantaneous model, a young age, however, they 
strongly overestimate the EWs of of Ha and H/3. Comparison of models cl and c2 shows 
that the adopted e-folding time for the host galaxy has a minor effect on tsF- 
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Table 4. Model results for each SF region 



Region 


Model 




tSF 


A/*,SF 


EW(Ha) 


EW(H/3) 


B - R 


AflB 




x 










(I0 b M ) 


s 

A 


A 


(mag) 


(mag) 




1 

J- 


inst 


10 


Myr 




20 


U 


89 

17. OZj 


f\9 

\J . UZj 


0.044 






2.5 Gyr 


1.46 


23 


7 


1 05 


1 7 

U.l 1 


389 




c 1 


7 


Myr 




996 

Zi C*U 


78 


84 

17. 


74 


029 

17. UZjC/ 




c2 


10 


Myr 


0.09 


237 


61 


0.82 


0.75 


0.032 




obs 








19.7 


5.4 








2 


inst 


8 


Myr 


n 1 7 


98 


7 


80 

17. OU 


60 


047 




rfl 


2.5 Gyr 


9 8 

Z* . O 


34 


1 


96 

17 . C/U 


0.42 


218 

17. Zj J. o 




r1 

k, _L 


7 


Myr 


1 1 


292 


77 


84 

17. 


73 

17. ( o 


024 

17. UZj^± 




c2 


10 


Myr 


0.09 


235 


60 


0.82 


0.74 


0.025 




obs 








28.9 


7.9 








3 


inst 


44 


Myr 


37.4 


7 


2 


0.60 


1.48 


0.227 




cO 






















2 


Gyr 


4Q 1 

4:i7 . _L 


1 08 


24 


56 
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- Col. 1: 


Star forming 


; region ID; col. 2: model; cols. 3, 


4: age and 


mass 


of the young stellar 



component; cols 5, 6: predicted equivalent width of the Ha and H/3 emission lines; col. 7: difference between 
the £?-band surface brightness of the older stellar population and that of the total stellar population predicted 
by the fit; col. 8: reduced x 2 of the fit. For each knot, the last line, labeled "obs", shows the observed Ha 
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and H/3 equivalent widths. 
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In summary, both population and evolutionary spectral synthesis models argue indepen- 
dently against the idea that the SF ring in Mrk 409 is forming continuously with a constant 
SFR since several 10 8 yr, as expected from the RR scenario. Instead, they consistently point 
to a recent (~10 Myr) enhancement of the SFR in all circumnuclear SF regions. 

Clearly, regardless of the model used, any conclusions related to the SFH of the ring 
should be considered with caution. It is known that both spectral synthesis approaches 
are plagued by ambiguities, especially when applied to complex emitting sources, such as 
starburst galaxies, involving multiple stellar populations of different age and metallicity, and 
being likely subject to different amounts of intrinsic extinction. 

Solutions by starlight are guided by the minimal linear combination of age vectors, 
i.e. they discard out of the available library SSPs with a low significance. Aside from the 
degeneracy problem, the solutions may therefore provide a very fragmentary view of the true 
SFH in a galaxy. In our concrete case, the absence of SSPs in the age interval between ~ 40 
Myr and ~ 1 Gyr should not be taken as foolproof of cessation of SF activities during that 
period, and the hypothesis of a low-level star formation in the ring cannot be rejected with 
certainty. 

Evolutionary synthesis models suffer as well from ambigui ties stemming not o nly from 



the age- metallicity degeneracy but, as thoroughly discussed in iGuseva et al.l (1200 If ) . by un- 
certainties in the intrinsic extinction and in the SFH assumed. The instantaneous and 
continuous model likely represent two extreme cases and rather poor approximations to the 
true SFH along the circumnuclear ring. For example, a continuous SF process over the last 
~ 1 Gyr with a varying or recently elevated SFR may also match the observations. How- 
ever, an investigation of more complex SFHs on the basis of the present data would likely 
not completely lift uncertainties in the interpretation of the SF ring. 

We stress that the satisfactorily agreement between the two conceptually different mod- 
els used is reassuring. We are therefore confident in that an instantaneous burst is a better 
description of the SFH of the individual SF knots in the ring. 

Furthermore, the morphology of the galaxy does not show any bar or non-a xisymmetric 



featur es, conditions required for the formation of a galactic RR. In a recent paper, IComeron et al 



(120091 ) note that Mrk 409 (NGC 3011) has no obvious bar, and place it into the group of 



galaxies that have "nuclear rings that cannot easily be explained in a resonance framework" . 

Hence, both the synthesis models results and the galaxy morphology speak against the 
presence of a RR in this galaxy, and favor the starburst-driven shock scenario as the most 
plausible to explain the origin of the SF ring in Mrk 409. 
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The nature of the nuclear region #3 is puzzling. As discussed already, its observed 
line ratios suggest an additional non-thermal source. Further support to this interpretation 
comes from spectral synthesis models which both fail to reproduce its observed SED without 
large systematic residuals (cf. Fig. [5]). Although inclusion of a power law component along 
with a variation in the extinction model might produce a better fit, we preferred not to do 
so, as this would probably not add solid constraints to the understanding of the nature and 
SFH of the nuclear region. 

However, assuming that the luminosity of region #3 is dominated by stars, we can 
formally check whether its energetic output is consistent with the blast wave scenario. Taking 
its predicted burst age and total stellar mass of t SF « 44 Myr and M+ |RP 3.7 x 10 7 M 



(Table HJ) at face value, we estimate from starburst99 flLeitherer et al.l Il999l . their Fig. 
Ill) the mean mechanical luminosity imparted to the ISM over tgp via radiative winds and 
SNe to ~ 3.8 x 10 41 erg s" 1 . This modest mechanical luminosity, which is to be taken as an 
upper limit, can indeed readily produce a super-shell of the observed dimensions within a 
few Myrs. 



From Eq. 2 of iMcCray fc Kafatosl (119871 ). the radius in pc of a shell writes as R s 



269(L3g/no) 1 / 5 ■ t^ 5 , where denotes the mechanical luminosity in units of 10 38 erg/s, 
t 7 = t/(10 7 yr) and no is the density of the ambient gas in cm . For L38 ~ 3800 the radius 
of the starburst-driven shell at the age tsF is R s ~3400/n ( 1 / 5 . Hence, the energetic output 
from the nuclear starburst region is sufficient for producing within t < tgp a super-shell with 
a R s between ~600 pc (inner SF ring) and ~2.3 kpc (outer SF ring), depending on n$. 

While the simple blast-wave scenario can formally be reconciled with the known prop- 
erties of Mrk 409, a major open question pertains to the actual SFH of region 7^3. 

While a conclusive assessment of this issue is out of reach with the present data, the 
co-evolution of a massive nuclear SF region and an AGN remains a viable and appealing 
interpretation for the overall characteristics of Mrk 409. 



5. Conclusions 



We presented PMAS integral field spectroscopy of the central 16" x 16" of the BCD 
galaxy Mrk 409. Mrk 409 is kno wn to undergo intense SF a ctivity within the central part 
of an old, elliptical host galaxy (jGil de Paz fc Madord 120051 ). Previous work (GdP03) had 
revealed a compact (diameter < 5") nuclear starburst region, surrounded by two nearly 
concentric rings of ionized gas emission with radii of 0.6 kpc and 2.8 kpc. The inner ring 
was proposed to be the result of in situ star formation in a dense shell of piled-up material 
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4000 4500 5000 5500 6000 8500 

wavelength (A) 

Fig. 9. — Observed and best-fitting synthetic spectrum with STARLIGHT (upper panel) 
and Pegase 2.0 (bottom panel) of the nuclear SF region #3. It can be seen from either dia- 
gram that a stellar population cannot fit the observed stellar continuum without systematic 
residuals. 
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on the surface of an expanding starburst-driven super-bubble. 

The main results from the present study may be summarized as follows: 

• The inner circumnuclear ring of Mrk 409 contains a young stellar population with a 
total mass of ~ 1.5 x 10 6 M & , which likely started forming almost coevally ~ 10 7 
yr ago. This is strong evidence to the interpretation that circumnuclear SF activities 
in Mrk 409 were triggered through the collision of a spherically expanding, starburst- 
driven super-bubble with the ambient denser material. Although the details of this 
process remain unclear, our study suggests that sequential star formation in a spherical 
shell, as observed in Mrk 409, may contribute considerably to the build-up of the young 
stellar component in BCDs. 

• While the circumnuclear SF ring presents properties typical of H II regions, diagnos- 
tic line ratios for the high-surface brightness starburst nucleus suggest an additional 
ionization source for this region, most probably a low-luminosity AGN. 

• The nuclear and circumnuclear SF components in Mrk 409 also differ considerably 
in their extinction properties: SF regions in the ring show a moderate and relatively 
uniform intrinsic extinction (0.2 < Cu/3 < 0.3), while the nuclear emission is subject to 
strong dust obscuration (Cu/3 ~ 0.9). Consequently, corrections for extinction based 
on the luminosity-weighted spectrum of the central region of a BCD can lead, in the 
case of Mrk 409, to a severe over-estimation of the Ha luminosity and star formation 
rate. 

• Mrk 409 shows a relatively high oxygen abundance (12 + log(0/H) ~ 8.40), without 
evidence for any significant chemical abundance gradient out to R ~ 0.6 kpc. 

• All the studied circumnuclear H II regions have electron densities in the low density 
limit. In the nuclear starburst region the electron density increases to a value as large 
as 500 cm" 3 . 

• The ionized gas component within the inner ~ 1 kpc of Mrk 409 displays a smooth 
kinematical pattern which is likely dominated by rotation; the total mass inside the 
SF ring is estimated to be ~ 1.4 x 10 9 M Q . 
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